Non-equilibrium polaron hopping transport through DNA 
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We study the electronic transport through short DNA chains with various sequences of base pairs 
between voltage-biased leads. The strong coupling of the charge carriers to local vibrations of the 
base pairs leads to the formation of polarons, and in the relevant temperature range the transport is 
accomplished by sequential polaron hopping. We calculate the rates for these processes, extending 
what is known as the P(i?)-theory of single-electron tunneling to the situation with site-specific 
local oscillators. The non-equilibrium charge rearrangement along the DNA leads to sequence- 
dependent current thresholds of the 'semi-conducting' current-voltage characteristics and, except 
for symmetric sequences, to rectifying behavior. The current is thermally activated with activation 
energy approaching for voltages above the threshold the bulk value (polaron shift or reorganization 
energy). Our results are consistent with some recent experiments.— 

PACS numbers: 71.38.-k, 05.60.-k, 87.14.gk, 72.20.Ee 



I. INTRODUCTION 

Experiments on long-range equilibrium charge transfer 
along DNA chains indicate that the dominant transport 
mechanism is polaron hopping of holes between the HO- 
MOs of adjacent base pairs in the DNA stack i^i^iii^il 
In most of these experiments a hole is injected into a 
guanine donor base by photo-chemical methods. The 
fraction of holes reaching other guanine base traps at 
a distance, separated from the donor guanine by vari- 
ous bridges is recorded. In the simplest case the bridges 
consist of a number of adenine bases, but more com- 
plicated bridges were investigated as well. The exper- 
iments showed a weak distance dependence for bridges 
longer than a few base pairs, which is in agreement with 
an activated hopping mechanism. Several theoretical 
articles argue that holes are localized on single (gua- 
nine) bases either by either solvation effects and/or struc- 
tural reorganizationi^ii^iii This localization can be in- 
terpreted as a polaron. The degree of localization is still a 
matter of debate^^, but many authors agree that confor- 
mational motion of the DNA is important to charge mi- 
gration in DNAfi^iiiii^ For short distances between gua- 
nine bases (number of intermediate adenine bases N < 3) 
the polaron migrates via superexchange (tunneling) be- 
tween the guanine bases, whereas for large distances the 
polaron undergoes hopping transport J^ii^iiSiii^ Several 
author a^°i^^'^^ have used an equilibrium polaron hopping 
picture to model results of various experiments. 

On the other hand, in non-equilibrium experiments, 
where short DNA chains were coupled to voltage-biased 
leads, different types of conduction were observed, reach- 
ing from quasimetallio21 via semiconducting2ii2^iS& to in- 
sulating behaviori^l This variance has not been explained 
(see, e.g., the review by Endres^S), but it suggests a 
strong influence of the environment and vibrations. Fur- 
thermore in several experiments^iSSi a strong tempera- 
ture dependence consistent with polaron transport was 



found (though another experiment"^*^ shows only weak 
temperature dependence). These observations suggest 
that polaron hopping plays an important role also in non- 
equilibrium transport through DNA. 

In most of earlier work dealing with electronic trans- 
port through DNA the vibrational excitations were ei- 
ther neglected altogetherj^ii^i^^ or they were included 
on a simplified level ^^ii^i^^'^ Other approaches describe 
polaron hopping as a classical biased random walk with 
hopping rates obtained by fits to experimen t "'^^'^^ or 
quantum chemistry calculationsj^^ sometimes including 
Coulomb charging effects,^" or they treat the problem nu- 
merically starting from the time-dependent Schrodinger 
equation4i In this paper, we study non-equilibrium 
polaron transport in short DNA strands connected to 
voltage-biased leads in the frame of rate equations for 
polaron hopping, with rates accounting for the transi- 
tions in the vibrational degrees of freedom as well. 

The DNA is modelled by a tight-binding chain where 
every site corresponds to either a guanine-cytosine (GC) 
or an adenine-thymine (AT) base-pair, which have dif- 
ferent on-site energies. Since both DNA strands have a 
direction (indicated by 5' and 3' at the ends) the local 
left-right symmetry is broken, and the hopping matrix 
elements depend not only on the combination of base 
pairs involved but also on the direction along the DNA 
molecule (compare table H]). Every base pair further cou- 
ples strongly to a local vibration, thus forming what is 
known as a 'small polaron'. Each vibrational degree of 
freedom, in turn, is coupled to an environment, which is 
responsible for dissipation and leads to a thermal occupa- 
tion of the vibrational states. We will evaluate the rates 
for polaron hopping in the spirit of what is known as the 
P(i?)-theory for electron tunneling in a dissipative envi- 
ronment modelled by a bath of oscillators42ii^'^ Here, 
instead of a bath of oscillators we have for each DNA 
base pair one localized vibrational mode which, however, 
is broadened due to the coupling to a dissipative envi- 
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ronment. 

Our main results are the following: (i) For DNAs 
with different sequences we observe pronounced differ- 
ences in the 'semiconducting' I-V characteristics, (ii) 
For (the usual case of) non-symmetric sequences we ob- 
serve strong rectifying behavior in the transport due 
to sequence-dependent hopping rates and occupation of 
the sites, (iii) For inhomogeneous sequences the current 
threshold positions are not directly connected to intrinsic 
energy scales, but depend on the charge rearrangement at 
finite bias, (iv) The current shows a strong temperature 
dependence following an Arrhenius law. The activation 
energy Ea is voltage dependent and approaches the bulk 
polaron value Ea ~ A/2 for high voltages (where A is 
the polaron shift or reorganization energy). 

II. MODEL AND TECHNIQUE 

We model the DNA with N base pairs by a minimal 
tight binding model, identifying every base pair with 
one electronic site. This is motivated by the fact that 
the molecular orbitals (HOMOs) of the charge carriers 
(holes) are located mainly on the purine bases, G or 
We describe the polaron hopping making use of 
the small-polaron theory, which has been developed to 
describe strong electron vibration coupling. The impor- 
tant point is that every DNA base can vibrate (quasi) 
independently from its neighbors, i.e. every site is con- 
nected to an independent oscillator. The vibrations in 
turn are coupled to a dissipative environment which al- 
lows for energy dissipation and relaxation to a thermal 
occupation. This environment can be modelled by cou- 
pling each oscillator to its own bath. 

Thus the Hamiltonian is H = H^i + i?L + ^^r + ^^tx + 

Ht,R + Hvih + -ffcl-vib + -H^bath, with 

i <ij> 
Ht,L + Ht,R = Y'^lr'^i + i"alc^r 

i ^ ^ 

i^oi-vib = ^ A, aja, (b, + B|) . (1) 

i 

The first term i/ci describes the electrons in the DNA 
chain with operators aj and ai in a single-orbital tight- 
binding representation with on-site energies of the 
base pairs and hopping between neighboring base 
pairs. Both on-site energies and hopping depend on the 
base pair sequence, e.g., the on-site energies differ for 
Guanine-Cytosine and Adenine-Thymine base pairs. For 
the direction-dependent hopping matrix elements we 
use the values obtained by Siebbeles et ai4^ who stud- 
ied intra- and inter-strand hopping between the bases 
in DNA-dimers by density functional theory. Adapting 



5'-XY-3'(all in eV) 
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0.119 


0.046 


-0.186 


-0.048 


C 


-0.075 


0.119 


-0.037 


-0.013 


A 


-0.013 


-0.048 


-0.038 


0.122 


T 


-0.037 


-0.186 


0.148 


-0.038 



TABLE I: Hopping integrals tij taken from Ref. |4^ and 
adapted to our model. The notation 5'-XY-3' indicates the 
direction along the DNA strand (see, e.g.. Fig. lb in Ref. [28|.) 



these results to our model of base pairs we obtain the 
hopping elements listed in tabled where, e.g., the num- 
ber in the row G and the column A denotes the hopping 
matrix element from a GC base pair to an AT base pair 
to its 'right', i.e., in the 3' directional 

The terms iJL,R(iiot written explicitly) with r = L, R 
refer to the left and right electrodes. They are modeled 
by non- interacting electrons, with operators cj^j. and Cj/r, 
with a flat density of states pe (wide band limit). The 
details of the coupling between the DNA and the elec- 
trodes are not the focus of this work. For our purposes it 
is sufficiently described by i?T,L + ^^t,r with tunneling 
amplitudes assumed to be independent of the base pair i 
and the quantum numbers of the electrode states v. The 
coupling strength is then characterized by the parameter 

The vibrational degrees of freedom of base i are de- 
scribed by i/vib, with bosonic operators e\ and Bi for the 
mode with frequency w^. The coupling of the electrons 
on the DNA to the vibrational modes is described by 
-f^ci-vib, where is the local electron- vibration coupling 
strength. Here we consider only the so-called stretch 
modes with frequencies hiOi — 16meV for a GC base pair 
and hiOi = llmeV for an AT base pair, which as shown 
by Starikov couple strongly to the electronsJ^ The cou- 
pling strengths are chosen in such a way that the reor- 
ganization energy or polaron shifts, Aa = 0.18 eV and 
Ag = 0.47 eV, fit the values extracted from experiments 
and listed by Olofsson et alM-. These values probably un- 
derestimate the effect of the solvent on the reorganization 
energy. 

The vibration of each base pair i is coupled to the local 
environment, i/i^bath, the microscopic details of which 
do not matter. It changes the vibrations' spectra from 
discrete modes tOi to continuous spectra, 

A(c^) y"rfte"^*0(t)({B|(t) + B,(i),B|+B,:}) 
f Vi{^) Vti^) \ 

(2) 

with frequency dependent broadening 77^(0;). The ac- 
tual form of rii{u!) depends on the properties of the bath. 
A reasonable choice which assures also convergence at 
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low and high frequencies is iji^oj) =^ 770 ^ ^(^^c — 
with r/o = 0.5 meV and a cutoff of the order of fkUc — 
0.045 meV. The couphng to the bath introduces the dis- 
sipation, which is crucial for the stability of the DNA 
molecule in current carrying situations where substantial 
amount of heat can be produced in the DNA. 

In order to describe the system with strong electron- 
vibration coupling we first apply the so-called polaron or 
Lang-Firsov unitary transformation 



H 



with the generator 



S 



- B] 



(3) 



(4) 



We introduce transformed electron and vibrational oper- 
ators, 



B, 



a-iXi 



Bi - - — a- a 



t. 



and polaron operators 



exp 



(5) 



Operators Xi with different indices i act on different vi- 
brational states, therefore they commute for all times. In 
terms of these quantities the Hamiltonian reads 



H =Ho + H' 

Ho = ^(ei - Ai)alai + ^ hoJi (bJb, + i 

i i ^ 



H' 



A,; = 



<ij> 



(6) 

(7) 

(8) 
(9) 



After these transformations we can proceed studying the 
effect of strong electron-vibration coupling in perturba- 
tion theory in H' . The small parameters are tij/Ai 
and V I which allows truncating the perturbation ex- 
pansion at lowest non-vanishing order in these parame- 
ters. From here on we will use the shifted on-site energy 
Q = — Ai in all expressions. 



Rate equation and current 

The small-polaron theory covers two limits of trans- 
port. At sufficiently low temperatures polarons form 



bands with bandwidth W 



where 



Wo exp 

Wq denotes the electronic bandwidth without vibra- 
tions.*^ At high temperatures the bandwidth W de- 
creases exponentially as the increasing number of multi- 
phonon processes destroy the coherence, and the band 
picture ceases to be valid. Transport is then accom- 
plished by a sequence of incoherent polaron hops. A 
rough estimate for the cross-over temperature is fc^T ~ 
huj [4 In {\/hw)]~^ ^ For the electron- vibration coupling 
strengths of interest in the present problem, room tem- 
perature is already well above this limit. 

To describe room-temperature transport it is there- 
fore sufficient to consider a rate equation for the diag- 
onal elements of the single particle density matrix, i.e. 

the occupation numbers of the sites pi{t) = (a\{t)ai{t) 

These occupation numbers evolve according to a master 
equation with transition rates which we obtain in an ex- 
pansion in H' from Fermi's Golden Rule. If we consider 
the rate for a hopping process from base pair (site) I to 
m, we have to take into account that also the vibrational 
states may change. If the initial and final states of the 
coupled system are denoted by / and F, the rates are 

Wz™ = f |(F \alxlaixi\l)\^ S{Ej - Ef).{10) 

In the following the vibrational states are not explicitly 
considered. Therefore, we trace out the vibrational de- 
grees of freedom Xi^m by summing over all initial vibra- 
tional states weighted by the appropriate thermal prob- 
ability and over all final state. Thus the transition rate 
from a state with site I initially occupied and site m ini- 
tially empty becomes 

Wi^^^\U„,\' J2 gi{xO\{xl\xi\x,)f 
x,,x: 



X - im + Ex, - Ex; + Ex„, - Ex^ ) , 

where Qi{Xi) is the probability of finding vibration I in 
state Xi . Rewriting the energy conserving delta- function 
by its Fourier transform we obtain 



1 



\tl„^? / dtei^''-'-'^'Pi{t)P^(t), (11) 



where 



Xi 



x. 



(12) 



=lCi exp 




dujDi{uj) 



dujDi {uj) 



\i \ cos (uj [t + ih(}/2]) 



sinh (fiw/3/2) 



coth(fiw/3/2) \ . (13) 
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The function Pi{t) is known from the '^P{E) theory", 
which describes tunneUng in a dissipative electromag- 
netic environment, modelled by an infinite set of oscil- 
lators. Here, instead of such a bath we have broadened 
local vibrational modes of two DNA base pairs involved 
in the hopping process. 

The calculation for the tunneling transition between 
the left (L) and right (i?) electrodes and the first or last 
site of the DNA chain ^ = 1 or Z = TV proceeds similarly, 
except that one has to trace also over the electrodes' 
electronic states, while we have to consider only the local 
vibration of the one site involved. Hence we have for the 
rates on the left, out and onto the DNA chain 



dE 
27rh 



(14) 



where V'^ — 27r|i'"ppe, fh{E) is the Fermi function in left 
lead, and Pi{E) is the Fourier transform of Pi{t). For 
the right interface a similar expression holds involving 
ME) and Pn{E). 

The master equation for sites in the DNA chain thus 
reads 

■^Pl = ^[ - Pli^ - Pin) yVlrn + (1 " Pi) PrnyVml , 
m 

(15) 

where the sum over m is restricted to nearest neighbors 
of For the base pair at the left end of the chain we 
get 

^^p,^-p,W^ + {l-p,)W]l 

+ [ - Pi (1 - P2) Wl2 + (1 - Pi) P2W21] , (16) 

and similar for the right interface. 

We are interested in the steady state, dpi/dt = 0, which 
develops for a constant applied bias. After solving the 
resulting self-consistent equations iteratively we can cal- 
culate the non-equilibrium current through the left lead. 



PiW^ + (1 - pi) wl 



(17) 



or for the right lead, which is the same since the current 
is conserved, /l — —Ik- 



where e™ and q are the on-site energies of base pairs m 
and Z, respectively. 

For large times, Pi{t) approaches a constant, 
lim Pi{t) = fCi. Therefore it can be separated into two 

t — >oo 

terms, one decaying in time and one constant: 

Pi{t)^Pi{t)+JCi . (19) 
Accordingly we can write 

Plit)Pmit) = Pl{t)P^{t) + ICMt) + IClPmit) + JCmlCl . 

The product Pi{t)Pm{t) describe transitions, where the 
number of vibrations changes on both sites, the next two 
terms describe changes in one of the two sites only, while 
the last term describes transitions without changes in the 
vibration state. When performing the time integration 
in eq. pT|) , this last term leads to a divergence when the 
two site energies are degenerate 



em) , (20) 



since in this situation the phenomena of resonant tunnel- 
ing occurs. In this situation the perturbation theory lim- 
ited to second order is not sufficient. Rather, one should 
sum up in a 'ladder'-approximation an infinite series of 
such terms, leading to a result with finite rates . ^^'^"^i^^ 

Alternatively, we can phenomenologically regularize 
the divergence of eq. (|20p by formally introducing an 
imaginary part to the level energies ei . This is motivated 
by the fact that they aquire a finite width due to the 
interaction with the vibrations or leads. In this way the 
hopping rates become finite. We further note, that the 
contribution of the vibration-free transitions (the con- 
stant term of eq. I20|) is multiplied by the factor ICiiCm- 
In our case, this factor is exponentially small. This cor- 
responds to the fact that we consider the limit where 
polaron hopping by far dominates polaron band trans- 
port. We therefore can ignore the terms oc JCmlCi in our 
analysis all together, i.e. we subtract them in eq. (|lip . 
The regularized hopping rates are therefore 



\tl,r. 



Pi{t)P^{t) - JCilC 



(21) 



Discussion of the hopping rates 



III. RESULTS 



For the hopping rates eq. fTTj) the situation differs from 
the usual P{E) theory: instead of one infinite vibrational 
bath each base pair (m and I) has its own vibration de- 
gree of freedom and we get products Pm{t)Pi{t), which 
become convolutions in energy space. The rates still sat- 
isfy detailed balance 



Wm; exp 



(18) 



Sequence effects 

Using the rate equation we will study now the charge 
transport and the non-equilibrium occupation of the sites 
for various DNA sequences. Both quantities depend 
strongly on the specific sequence. All DNA sequences 
are 'semi-conducting' since the Fermi energy lies in the 
HOMO-LUMO gap, i.e. weU above the HOMO states 
which carry the transport. 
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FIG. 1: (Color online) I-V characteristics for two DNA 
strands with sequences 5'-GGGGGGGG-3' (dash-dotted line) 
and 5'-GAAAAAAG-3' (solid hne) with the following pa- 
rameters: Base pair on-site energies ~ — 0.26eV, ec = 
-1-0.25 eV, polaron shifts Aa = 0.18 eV and Ac = 0.47 eV, 
Fermi energy _Ef ~ OeV, symmetric coupling to leads with 
linewidths Fl = Fr = 0.01 eV, vibrational energies huiA = 
llmeV, hiJG ~ 16meV, and room temperature ksT = 
25 meV. The inset shows the absolute value of the current 
on logarithmic scale. The current for the second sequence 
shows rectification by a factor of ~ 200. 



FIG. 2: (Color online) Differential conductance (logarithmic 
scale) as a function of applied bias for five different DNA se- 
quences with parameters as in Fig.[T] For the homogeneous se- 
quences the threshold, i.e. the position of the maximum of the 
difi'erential conductance, is set by on-site energy of the consid- 
ered base pairs. For the inhomogeneous sequences, however, 
the threshold is not determined by the internal energy scales 
alone. The sequence-dependent thresholds lie in between the 
limits set by the homogeneous sequences. For some sequences 
(e.g. 5'-GAAAAAAG-3'), the peaks are broadened due to 
'renormalized' tunneling. 



Figure [T] shows the I-V characteristics for two such 
sequences, 5'-GGGGGGGG-3' (green, dash-dotted Hne) 
and 5'-GAAAAAAG-3' (black, sohd hne). The first se- 
quence displays the 'semi-conducting' behavior with a 
gap characterized by the distance of the Fermi energy to 
the onsite energy of the G base (shifted by Aq). Due to 
its electronic symmetry the I-V characteristic is symmet- 
ric with respect to the applied bias. On the other hand, 
the second DNA sequence shows strong rectifying behav- 
ior, despite of its seemingly symmetric sequence. The 
reason for this asymmetry lies in the electronic asymme- 
try of the hopping amplitudes, together with the incoher- 
ence of the hopping processes between DNA base pairs. 
This can easily be understood: For positive bias the hop- 
ping 'bottleneck' of the system is at the crossover from 
A to G at the 3' end of the strand. There, the polaron 
needs to overcome an energy barrier mediated by vibra- 
tional excitations. For negative bias the 'bottleneck' is 
at the crossover from A to G at the 5' end of the strand. 
Due to the opposite direction of the dominating hopping 
process, with I^gaI > |*ag| (compare Table HI, the cur- 
rent for negative bias is higher than for positive bias. 
Thus, inhomogeneous sequences will in general display a 
rectifying, semi-conducting I-V characteristic. The rec- 
tification effect will be weaker for longer and more dis- 
ordered sequences, as more 'bottlenecks' in either direc- 
tion appear. Note that no rectifying behavior would be 
observed if we model the transport as a coherent tran- 
sition through the total length of the chain ( 'Landauer 
approach ' ) Mi^^^ 



We now study the sequence dependence of the current 
threshold, or equivalently the postion of the peak in the 
differential conductance dl/dVh (both differ only by a 
term proportional to temperature) . Figure [H shows the 
differential conductance as a function of the applied bias 
for 5 different DNA sequences 5'-AAAAAAAA-3', 5'- 
GAAAAAAG-3', 5'-GGAAAAGG-3', 5'-GGGAAGGG- 
3', and 5'-GGGGGGGG-3'. For the homogeneous se- 
quences the threshold is equal to the on-site energy of the 
considered base pairs (el4 = 2eA for 5'-AAAAAAAA-3' 
and eVb = 2?g for 5'-GGGGGGGG-3'). For the inhomo- 
geneous sequences the threshold lies in between the limits 
set by the homogeneous sequences, i.e. it is not deter- 
mined by the internal energy scales alone. The varying 
threshold is a consequence of the way the charges are 
rearranged along the DNA molecule, which of course is 
very sensitive to the considered sequence. 

As discussed above, eqs. (|14p describe the tunneling 
rate from the electrode to the adjacent DNA base pair. 
The tunneling is modified by the vibrational modes which 
can be excited, depending on the applied bias. This 
'renormalized' tunneling can lead to a very broad differ- 
ential conductance peak which is very different from the 
usual (derivative of) Fermi function form, as observed 
most prominently for the sequence 5'-GAAAAAAG-3'. 
Note that there is nearly no modification on the low bias 
side of the peak. 
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FIG. 3: (Color online) I-V curves for the two sequences 5'- 
GAAAAAAG-3' (solid line) and 5'-GGAAAAGG-3' (dashed 
Une) (parameters see Fig. [1]). Despite the very similar se- 
quences, the I-V show clear diflterences. The inset shows the 
local chemical potential "1> for the last guanine base (at the 3' 
end) for both sequences at various bias voltages. Equivalent 
behavior between local potential and I-V is visible 
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FIG. 4: (Color online) Local chemical potential $ for all base 
pairs of the DNA strand with sequence 5'-GAAAAAAG-3' 
(black lines) and 5'-GGAAAAGG-3' (red lines) at various bias 
voltages and with parameters as in Fig. [T] The local poten- 
tial drops differently for the two sequences, implying different 
conduction properties. 



Local chemical potential 

As discussed above, the I-V characteristic of a DNA 
molecule is affected by bias and sequence dependent 
charge rearrangements on the DNA base pairs. For the 
ease of displaying these effects, i.e, both small deviations 
from an occupation 1, as well as occupations near 0, we 
introduce a local chemical potential <i>i, defined by 



fcRTln 



1 



1 



(22) 



This quantity is superior to the occupation in visualizing 
the non-equilibrium charge rearrangement, because it re- 
acts sensitively to even small changes in the occupation. 

Figure [3] shows the I-V curves for the two DNA 
molecules 5'-GAAAAAAG-3' (black, sofid fine) and 5'- 
GGAAAAGG-3' (red, dashed line), and the inset shows 
the local chemical potential $ for the last guanine base 
(at the 3' end) for both sequences. Although the se- 
quences are very similar, the I-V characteristics differ 
strongly in the maximum current and in the way the 
current increases for increasing bias voltage. The current 
of the second sequence has reached a plateau already at 
about Vh = —0.8 V, whereas the black curve has not lev- 
eled off even for Vb = —1.5V. This strong deviation from 
a Fermi function behavior is in part a consequence of the 
renormalization of the tunneling rates by the vibrations. 

This difference in the I-V characteristics is reflected 
in the local chemical potential most prominently at 
the last guanine base of both sequences, as shown in the 
inset. At low bias both sequences behave in the same 
way: the potential increases equally with the applied 
bias. The DNA is not conducting and therefore, the sit- 
uation is similar to the charging of a capacitor. At the 



drop-off around Vb = —0.3 V, the current sets in and 
a potential drop between base pair and lead is estab- 
lished. In correspondence to the current, the local chem- 
ical potential for the second sequence 5'-GGAAAAGG-3' 
levels off, whereas the potential of the first sequence 5'- 
GAAAAAAG-3' never reaches a plateau in the range up 
to Vb = -1.5 V. 

To give a feeling for the total charge rearrangement 
Figure [3] shows the local chemical potential $i of the 
two DNA sequences for all base pairs i and all voltages. 
The chemical potential landscape also suggests how the 
bias voltage Vb applied to the leads drops over the en- 
tire DNA molecule. Regions of good conductivity show 
almost no voltage drop, as seen for the stretches of ade- 
nine bases in the middle of both sequences. On the other 
hand most of the voltages drops at the base pairs close 
to the interfaces. The potential/voltage drop over for 
the entire sequence 5'-GAAAAAAG-3' is less than for 
5'-GGAAAAGG-3'. This suggests that the the first se- 
quence is better conducting than the latter, which is in 
accordance with their I-V curves (Fig. [3]). 



Temperature dependence and activation energy 

In the experiments of Ref. [l] the current through 
bundles of long homogeneous DNA molecules showed a 
strong temperature dependence. The data could be fit- 

with 



ted by an activation law liV) = a{V) exp 

a voltage dependent prefactor a(V) that also shows a 
temperature dependence for the case poly(dG)-poly(dC) 
bundles. Asai^^ has used the Kubo formula for a polaron 
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FIG. 5: (Color online) Activation energy Ea for polaron hop- 
ping of a homogeneous DNA with 15 G-C base pairs as a 
function of the polaron shift A for voltages Vb ~ 0.04 V (solid 
line), Vb = 0.55 V (dashed line), and Vb = 0.8 V (dash-dotted 
line). All other parameters as in Fig. [l] For comparison, the 
dotted line shows the activation energy of polaron hopping 
conduction in bulk at high temperatures {Ea ~ 0.5A).^^ 



hopping model to obtain a similar relation for the linear 
response conductivity. 

Our results are obtained in a non-equilibrium situation 
and also show a strong temperature dependence. An Ar- 
rhenius plot of the current vs. temperature shows linear 
behavior, indicating that the current is indeed an acti- 
vated quantity (though we also observe deviations from 
a perfect Arrhenius law). Fitting the temperature depen- 
dence of our data by an Arrhenius law allows us to esti- 
mate the activation energy for a given bias voltage and 
polaron shift A. Figure [5] shows the activation energy Ea 
obtained by this fitting as a function of the polaron shift 
at three different bias voltages for a homogeneous DNA 
strand with 15 GC base pairs.—- The activation energy 
is proportional to Aq, but the proportionality factor dif- 



fers depending on the applied bias voltages. For voltages 
smaller than the gap, the activation energy also includes 
the energy needed to overcome the gap. For voltages 
above the threshold the proportionality factor between 
activation energy and polaron shift is about 1/2, consis- 
tent with the high-temperature value for the bulk polaron 
hopping conduction Ea = 0.5A (see green/dotted line in 
Fig.li])^ 



IV. SUMMARY 

We have investigated the non-equilibrium polaron hop- 
ping transport in short DNA chains with various se- 
quences, coupled to voltage-biased leads in the frame of 
rate equations which take into account inelastic transi- 
tions in the local vibration degrees of freedom. Our the- 
ory is formally an extension of the so-called P{E) the- 
ory of tunneling in a dissipative electromagnetic environ- 
ment. We find semi-conducting I-V characteristics with 
thresholds that are very sensitive to the considered DNA 
sequence. For all non-symmetric sequences (which is the 
typical case) we observe rectifying behavior (Fig.[T]). The 
sequence dependent thresholds are not directly connected 
to intrinsic energy scales (Fig. [5]), rather they are inti- 
mately related to the non-trivial charge rearrangement 
along the DNA molecule at finite bias. We have visual- 
ized this effect by displaying the local chemical potential 
$i (Fig. 0]). As expected for polaron hopping, the cur- 
rent is thermally activated with a temperature depen- 
dence following an Arrhenius-law. The activation energy 
Ea is voltage dependent and approaches the bulk polaron 
value Ea — 1/2 A (A: polaron shift) for voltages above 
the threshold (Fig. [S]). 
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